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■ Abstract 

\o 

I Semiclassical approximations often involve the use of stationary phase approximations. This 

method can be applied when Ti is small in comparison to relevant actions or action differences in 
the corresponding classical system. In many situations, however, action differences can be arbitrar- 
ily small and then uniform approximations are more appropriate. In the present paper we examine 
different uniform approximations for describing the spectra of integrable systems and systems with 
mixed phase space. This is done on the example of two billiard systems, an elliptical billiard and 
a deformation of it, an oval billiard. We derive a trace formula for the ellipse which involves a uni- 
form approximation for the Maslov phases near the separatrix, and a uniform approximation for tori 
of periodic orbits close to a bifurcation. We then examine how the trace formula is modified when 
the ellipse is deformed into an oval. This involves uniform approximations for the break-up of tori 
and uniform approximations for bifurcations of periodic orbits. Relations between different uniform 
approximations are discussed. 
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1 Introduction 



Semiclassical trace formulas are important tools for analyzing spectra of quantum systems. They have 
found a wide range of applications in recent years (see for example |l|, ^, Q[). Trace formulas have 
been derived for integrable and chaotic systems when periodic orbits lie on tori in phase space or are 
isolated M,fh ||, 0, Families of orbits in systems with more general symmetries have also been 
treated || [l^] . 

Most systems, however, are neither integrable nor chaotic, but have a phase space which is in- 
tricately divided into regular and chaotic regions. In these systems there exist classical structures 
on all scales, and neighbouring periodic orbits can have arbitrarily small action differences. For the 
semiclassical approximation this has the consequence that in many cases stationary phase approxima- 
tions cannot be applied since they rely on the fact that action differences of periodic orbits are large 
in comparison to fi. Instead one has to use approximations which are uniformly valid in two small 
parameters, Ti and e, where e describes the separation of neighbouring periodic orbits. 

A typical situation where small action differences of periodic orbits occur is near bifurcations of 
periodic orbits. Bifurcations are a ubiquitous phenomenon in mixed systems. They occur any time the 
stability angle of a stable orbit is a rational multiple of 2ir. If, for example, one changes an external 
parameter of a system by an arbitrarily small but finite amount (or the energy in a generic system), 
then in general an infinite number of bifurcations occur, most of them for long periodic orbits. 

Bifurcations occur in different forms, but the number of generic forms in two-dimensional systems 
is finite. They were classified by Meyer 11] and Bruno [l3|]. Their form depends on the repetition 
number m of an orbit that bifurcates. The semiclassical treatment of these generic bifurcations 
was investigated by Ozorio de Almeida and Hannay 14]. They derived contributions to the trace 
formula from periodic orbits near bifurcations in terms of standard diffraction catastrophe integrals. 
These approximations have for example been applied for treating tangent bifurcations and pitchfork 
bifurcations ||T5| , [T^| . For the largest class of bifurcations with m > 4 the results of Ozorio de Almeida 
and Hannay were extended in [17] by including higher order terms in the normal form expansion for 
the bifurcation. In this way a uniform approximation was obtained which interpolates over the regime 
from the bifurcation up to regions where Gutzwiller's approximation can be applied. 

Another situation in which small action differences occur is in the quasi-integrable regime, i. e. for 
small perturbations of integrable systems. Due to the perturbation all rational tori of the integrable 
system break up into a finite number of periodic orbits which have small action differences when 
the perturbation is small. Ozorio de Almeida derived a uniform approximation which is valid if the 
splitting of the orbits is small [jD| . A formula for the break-up of families of orbits due to more general 
symmetries was derived by Creagh [19]. Tomsovic, Grinberg and Ullmo extended the result of Ozorio 
de Almeida and obtained a uniform approximation which interpolates between the torus contribution 
and Gutzwiller's approximation |2(], 21]. 

Small action differences of neighbouring periodic orbits can occur also in integrable systems, since 
bifurcations of periodic orbits can occur also there. These bifurcations do not belong to the class 
of generic bifurcations, since they result in the appearance of whole new tori of periodic orbits. A 
uniform approximation for the contributions of tori that result from the bifurcation of a stable orbit 
was derived by Berry and Tabor ||] and Richens [ 22 ] . 

The present paper contains an investigation and discussion of the various concepts which are 
mentioned above. This is done on the example of two billiard systems, an ellipse and a deformation 
of it, an oval. The classical motion in an ellipse is integrable and has several characteristic properties. 
There is a separatrix which separates two kinds of motion, around the two foci and between the two 
foci. Furthermore, in addition to the usual tori of periodic orbits there are also two isolated periodic 
orbits in the ellipse (and their repetitions), one stable and one unstable. When the eccentricity of the 
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ellipse is increased, then new tori of periodic orbits arise through bifurcations of the stable orbit and 
its repetitions. All these classical properties have an influence on the semiclassical approximation. 
After a brief review of the classical and quantum mechanics in the ellipse we discuss the semiclassical 
EBK-quantization for the ellipse. Near the separatrix this quantization condition has to be modified 
by a (non-integer) uniform approximation for the Maslov index, since this index is discontinuous at 
the separatrix which leads to ambiguities and inaccuracies in the semiclassical quantization p4| ]. 
From this modified EBK-quantization condition we derive a trace formula for the ellipse. We show in 
detail, how the contributions of regular tori of orbits, of tori near bifurcations, of the isolated periodic 
orbits and the area and perimeter contributions can be obtained from the modified EBK-quantization 
conditions. We close this section with a numerical examination of the trace formula. 

We then examine an oval-shaped billiard which can be considered as a deformation of an ellipse. We 
investigate how the semiclassical approximation for the ellipse is changed when the ellipse is deformed 
into the oval. Due to this perturbation all tori of periodic orbits break up into a finite number of 
orbits. When the tori are not close to a bifurcation then the uniform approximation for the break-up 
of tori can be applied. For tori near bifurcations one has to use a different approximation. It will be 
shown that the (slightly modified) uniform approximations for generic bifurcations describe also the 
break-up of the tori near bifurcations. We present a detailed numerical examination of semiclassical 
approximations in the oval, and we discuss the relations between different uniform approximations 
that are used in this paper. 



2 The elliptical billiard 

We discuss in this section various properties of the elliptical billiard. This is an integrable billiard 
system whose boundary is defined by 

i + v i = 1 ' « 

where a and b < a are the lengths of the semi-major and semi-minor axis, respectively. The ellipse 
has two focal points with coordinates (±c, 0), where c = V a 2 — b 2 , and its eccentricity is defined by 
e = c/a. The billiard area is A = nab and its perimeter is L = 4aE(e) where E(k) is the complete 
elliptic integral of the second kind |25|] (see Eq. (|ll]) below). 

The classical properties in the elliptical billiard have been investigated under different view points. 
This includes the treatment of the action-angle variables p6fl, the caustics of the classical motion 



[pq , 27(1 , Poncelet's theorem |2g, 2J|, the billiard map []3C| plf] , and the periodic orbits [29| |3^, |33[. A 
treatment of the Schrodinger equation for the elliptical billiard can be found in fl34| . In the following 
we briefly review classical and quantum properties of the elliptical billiard because they are used for 
the semiclassical approximation. We then discuss the EBK-quantization for the billiard and a uniform 
extension of it, and we derive a trace formula in terms of the periodic orbits of the system. 

2.1 Classical mechanics of the elliptical billiard 

The classical motion of a particle in an elliptical billiard is conveniently described in elliptical coordi- 
nates 

x = c cosh u cos v 

y = c sinhu sinv , (2) 
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where u is restricted to < u < arctanh(6/a) and v is a periodic coordinate with period 2ir. In terms 
of these coordinates the Lagrangian has the form (with mass m = 1/2) 



the canonical momenta are given by 



Pu 



Pv 



— (cosh 2 u - cos 2 v )(n 2 + v 2 ) , 



-— = —it (cosn u — cos v) 

CJ it 

dL ° 2 ■ ( V.2 2 x 

-z— = — w (cosh u — cos v) , 

CJ V 



(3) 



and the Hamiltonian is 



H 



pI+pI 



c 2 (cosh u — cos 2 v) 



(4) 



(5) 



The two conserved quantities of the system are the energy and the product L1L2 of the two angular 
momenta about the two foci 

p 2 sinh 2 u — pi sin 2 v 



1-^2 



(6) 

cosh u — cos^ v 

It is more convenient to use instead of L1-L2 another conserved quantity which is energy independent 
and is determined by the geometrical properties of a trajectory only. This is a = L1L2/E whose 
values are restricted to the range b 2 > a > —c 2 . The upper limit a = b 2 corresponds to the motion 
along the boundary and the lower limit a = —c 2 corresponds to the motion along the minor axis. 
In terms of E and a the canonical momenta are given by 



Pi 

pI 



= E(c 2 sinh 2 u — a) 
= E(c 2 sin 2 v + a) . 



(7) 



There are two different kinds of motion in the ellipse depending on the sign of a. For < a < b 2 the 
trajectories have a caustic in form of a confocal ellipse with semi-minor axis b' = \fa. The motion goes 
around the two focal points and is composed of a libration in the coordinate u and a rotation in the 
coordinate v. For — c 2 < a < the caustic of the motion is a confocal hyperbola with semi-conjugate 
axis b' = yj—cx. and semi-transverse axis a' = V c 2 + a. The motion is composed of a libration in the 
coordinate u and a libration in the coordinate v, and the trajectories cross the x-axis always between 
the two focal points. Both kinds of motions are separated by a separatrix which consists of orbits with 
a = that go through the focal points of the ellipse. 

The periodic orbits of the elliptical billiard can be determined by introducing action-angle variables. 
The actions are given by 



— (fp u du = I du J c 2 sinh 2 u — a 
2ir J 7r J UQ 

1 r , 2^/e r>i / 



where 



c 2 sin 2 v + a , 



for a > : no = arcsinh— , u\ = arcsinh- , Vq = , v\ = — , 

c c 2 

b b' it 

for a < : no = , u\ = arcsinh- , no = arcsm — , v\ = — , 

c c 2 



(8) 



(9) 
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and b 1 = \J\a\. We choose here a definition of I u which is actually the action for only half a cycle if 
a < 0. The reason is that otherwise I u is discontinuous if a changes sign. We remark that the actions 
in a half-ellipse which is cut along the x-axis have the values given by Eqs. (|8|) and (||). 
The integrations in (||) yield 
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for a > 
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(10) 



where k = c/a' = c/\/c 2 + a, and the functions 

f da 

F(^«) 



\/l — k 2 sin 2 a 



E(</?, k) = f \fl — k 2 sin 2 a da , 
Jo 



(11) 



are the elliptic integrals of the first and second kind [|25| . The quantity k is called the modulus of the 
elliptic integrals. The corresponding complete elliptic integrals are denoted by E(k) = E(7r/2, k) and 
K(k) = F(tt/2,k). 

The tori of periodic orbits of the elliptical billiard are determined by the condition that the quotient 
of the angular frequencies u u and lo v is rational: 
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(12) 
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(13) 



one obtains the conditions for periodic orbits. We will consider in the following the cases a > and 
a < separately. 



The case a > 0: 

The condition for periodic orbits is 



F (arcsin 
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b 2 



2m ^,7T , 2m . . 
— F (-,« ) = — K(«) 
n z n 



(14) 
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This equation has a solution for all n = 3,4, . . . and 1 < m < n/2, and the integers n and m are the 



number of reflections of an orbit and its rotation number, respectively. Eq. (14) determines the value 
of a for a family of periodic orbits. 

We briefly discuss the case m = n/2. For this case there exist also periodic orbits which are the 
orbit along the major axis and its repetitions. However, these orbits are isolated and do not appear in 
families. This is reflected by the fact that for m = n/2, there exists no solution of Eq. ( |i~4| ) (although 
it does exist for Eq. fli~2j)). In the limit a — > 0, the first argument of the elliptic integral on the right- 
hand side of Eq. ( |14| ) approaches vr/2, but the integral itself is divergent in this limit since k — > 1. The 
correct limiting behaviour can be obtained from the relation [p5[ 



F(<p, k) + F(V>, k) = F(— , k) if \/l — k 2 tan (p tan ip = 1 . 



(15) 



From this relation follows that 
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it b 2 - 

F (2' K ) - F (arcsnW — — 
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F(arctan - , 1) = log 

b c 



(16) 



A different way of writing condition (fbS) is obtained by inverting the equation which leads to 



'b 2 — a 2m 
— = SU{ — K) ' 



(17) 



where sn(it) is a Jacobian elliptic function. In the above notation, the modulus k, which acts as an 
independent variable, is omitted. 

The lengths of the periodic orbits are given by 
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(18) 



where Eqs. (10) and (14) have been used and Z(u) is Jacobi's zeta function 

E(vt/2,k) 



Z(u) = E((/3, k) — F(<p, k) 



F(n/2, K ) 



(19) 



Here u = F((p, k) and the dependence of Z{u) on the modulus k is again omitted in the notation. 

Addition theorems for Jacobian elliptic functions and Jacobi's zeta function can be used to obtain 
algebraic expressions for lengths and a-values of periodic orbits. For example, for the case n = 4 and 
m = 1 one obtains 



sn — 

v 2 ; 



1 



i + Vi 



y 2> 21+ vT 



(20) 



From this follows that the periodic orbits have the conserved quantity a = b 4 /(a 2 + b 2 ) and length 
I = Wa 2 + b 2 . 



6 



The case a < 0: 

For a < the condition for periodic orbits is 



F(arcsin 



62 



b 2 



a k 



^F(|,I) 
n 2 k 



2m 1 

n k 



(21) 



with the restriction that n has to be an even integer, since the bounces at the boundary occur 
alternately in the upper and lower half of the ellipse. This requirement is a consequence of the 
definition of I u in Eqs. (g) and @ which is the action for only half a cycle if a < 0. 

The tori of periodic orbits can be labelled by n = 4, 6, . . ., and 1 < m < n/2, where the integers n 
and m are the number of reflections of an orbit and its libration number, respectively. Eq. ( ^T|) has a 
solution for all these values of n and m in the range < a < — oo. However, if the value of a for a 
solution is smaller than — c 2 , then the torus is complex, as is the modulus 1/k. If one decreases the 
ratio b/a then the torus becomes real. This happens when 



b mir 
— = sin 

a n 



(22) 



which follows from Eq. (21) with a = — c 2 and F(y?,0) = (p. For a = — c 2 the torus has zero extension 



and coincides with the orbit along the minor axis. Thus all tori arise through bifurcations of the stable 
orbit (and its repetitions) when the ratio b/a is decreased from its starting value b/a = 1 in a circle. 

We discuss again the case m = n/2 which corresponds to the periodic orbit along the major axis 
and its repetitions. Also for a < there exists no solution for Eq. (f2~l|) if m = n/2 as can be seen by 
using (JT5|) . This yields 
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The inverted form of condition d2l]) is given by 
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(24) 



where the modulus of the functions sn and K is now 1/k. 
The lengths of the periodic orbits are given by 
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(25) 



where Eqs. (10) and (pl|) have been used. 

We give again explicit expressions for the case n = 4 and m = 1. The values of sn(K/2) and 
Z(K/2) are now given by ( |20| ) with k replaced by 1/k. It follows that the periodic orbits have the 

-6 4 /(a 2 



conserved quantity a 



b ) and length I 



4a 2 /Va^ 



b 2 . 



Detailed illustrations of properties of periodic orbits in the ellipse can be found in [32|. 
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2.2 Quantum mechanics of the elliptical billiard 

The Schrodinger equation in elliptical coordinates has the following form (in dimensionless units 
% = 2m = 1) 

( d 2 d 2 \ 

— I + ^2 ] ^( U > -U ) = -^^(cosh 2 u - cos 2 v)t/)(u, v) . (26) 
Writing tp(u,v) = ipi(v)ip2(u) this equation separates into two ordinary differential equations 

ip"(v) + (b - h 2 cos 2 v)fa(v) = (27) 
-V4>) + (fr-^ 2 cosh 2 u)^ 2 (u) = 0, (28) 

where /i = \/l£c and b = E(c 2 + a). We use here and in the following the notation of Morse and 



Feshbach [34]. The first equation ( P7|) is Mathieu's equation. For every value of h there is a countable 
number of values of b for which it has periodic solutions of period tt or 2tt. These values are denoted by 
be r (/i), r = 0, 1, . . ., and ho r (h), r = 1,2,.... The corresponding solutions are the Mathieu functions 
Se r (h, z) and So r (h, z), where z = cosv, < v < ir. They are even and odd, respectively, with respect 
to reflection on the x-axis {v — * —v). 

The second equation (|2^) is Mathieu's equation for an imaginary argument: With the values be r 
and bo r for the constant b, the solutions of this equation that are regular at u = are the radial 
Mathieu functions of the first kind Je r (h, z) and Jo r (/i, z) where z = coshti. The energy eigenvalues 
of the Schrodinger equation follow from the condition 

Je r (h, cosh U) = 0, r = 0, 1, . . . 

Jo r (/i,cosh[/) = 0, r = l,2, ... (29) 

where U = arctanh(6/a). The functions Je and Jo have the properties [dJe r (h, cosh u)/du] u= o = 
and Jo r (h, cosh u)\ u =o = 0. 

The solutions can be separated into the four symmetry classes of the elliptical billiard. These 
symmetry classes will be denoted by two letters, where the first denotes the boundary condition on 
the x-axis and the second that on the y-axis. For example, DN denotes the symmetry class of wave 
functions which satisfies Dirichlet boundary conditions on the x-axis and Neumann boundary condition 
on the y-axis. 

The solutions of the Schrodinger equation for all four symmetry classes are, up to a normalization 
constant, 

DD : So2r+2(ft-i cos v) Jo2r+2(h, cosh u) DN : So2r+i (h, cos v) Jo2 r .fi (h, cosh u) , , 
ND : Se2r+i (h, cos v) Je2 r +i (h, cosh u) NN : Se2 r (h, cosv) Je2r(7i, cosh u) , 

where r = 0, 1, . . .. In the limit h — ► these solutions reduce to the solutions of the circular billiard, 



i.e. to a product of trigonometric and Bessel functions [34|. 

A numerical examination of the energy spectrum of the elliptical billiard in dependence on the 
ratio of the two half-axis can be found in l36|l . 



2.3 The semiclassical approximation for the elliptical billiard 

In this section we derive a semiclassical trace formula for the spectral density of the elliptical billiard. 
We follow the method of Berry and Tabor || and start with the EBK-quantization. The semiclassical 
spectral density which is obtained from it is then reexpressed by applying the Poisson summation 
formula to it. From this representation the trace formula for the ellipse is derived. 
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2.3.1 The EBK-quantization 



The EBK-quantization conditions for the elliptical billiard have been examined in [26]. They are given 
by 



a > : 
a < : 
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rit, = 0, 1, . . . 



(31) 



The condition for /„ and a > corresponds to periodic boundary conditions, and the factor of 2 
which appears in the condition for I u and a < is due to the fact that I u is the action for only 
half a cycle. The formulation of the EBK-conditions for the full ellipse is not very convenient. If the 
a- value of a semiclassical state changes sign (as a consequence of changing the ratio b/a), then the 
state is described by different quantum numbers than before. For that reason, it is advantageous to 
formulate the semiclassical quantization for a half-ellipse where this problem does not appear. The 
semiclassical levels of the full ellipse are then obtained by adding the spectra of two half-ellipses which 
have Dirichlet (D) or Neumann (JV) boundary conditions on the x-axis, respectively, (and Dirichlet 
boundary conditions on the remaining arc). For these two systems the EBK-conditions have the form 



(32) 



As one can see, the quantization conditions are discontinuous at a = 0, and as a consequence the 
semiclassical quantization is inaccurate near a = 0. This was examined in detail in [23, 24]. It 
was found that near a = there are sometimes two different semiclassical approximations for one 
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quantum state, and sometimes there is no semiclassical approximation. A remedy to this problem is 
the application of a uniform approximation near a = 0, which yields quantization conditions in terms 
of a noninteger Maslov index that interpolates smoothly between the cases a < and a > 0. This 

and 



uniform approximation has been derived in [23] by expanding the potential terms in Eqs. (f27j) and ( |2q ) 
up to second order in u and v. The solutions of the corresponding differential equations are parabolic 
cylinder functions. By matching the asymptotic form of these solutions with the EBK-solutions one 
obtains a uniform approximation for the Maslov phase. As a consequence the conditions ( |32[) are 
replaced by 



, 1 , fo(a) 2/? A (-q) 
D : I u = n u + --\ J v = n v -\ n u , n v = 1, 2, . . . 

2 7T 7T 
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1 T,/ 1 

-argr(- + ? a) 



(34) 



6 = arctanf Vl + e 2 ™ + e™) , (35) 

and (3a(&) is given by (34) with the replacement of 6 by tt — 8. The original conditions ( |32[) are 
recovered with the limiting values of (3s and /3a 



a -> +oo : (3 s ,(3a 



TT 
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-oo : ^5^0, /3a 



7T 
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(36) 
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With this uniform approximation for the Maslov phase the previous ambiguity of the semiclassical 
quantization is removed. A further advantage of the uniform approximation is that it removes a 
semiclassical degeneracies of energy levels with a > 0. According to the quantization condition ( |32| ) 
the even and odd states {D and N) have the same energy values (except for the lowest even state). 
This degeneracy is a semiclassical degeneracy since the true quantum levels are not degenerate. 

An alternative possibility for treating the semiclassical influence of the separatrix consists in the 
inclusion of complex tunneling orbits. This is done for the ellipse in fl37fl . 



2.3.2 The spectral density 



We derive now a trace formula for the spectral density of the elliptical billiard in terms of the periodic 
orbits of the system. We follow the method of Berry and Tabor Q and apply the Poisson summation 
formula to the semiclassical spectral density which is obtained from the EBK-quantization. We do 
this here for the two cases of a half-ellipse with Dirichlet or Neumann boundary conditions on the 
x-axis. The spectral density is given by 



d(E) 



(37) 



n u ,n v =0 



00 roo roc y y 

= z2 / dl u / dl v exp{2irin(I u - + 2mm(I v - -^)} 8(E - E n n J , 

where after the application of the Poisson summation formula the integration variables have been 
changed from n u and n v to I u = n u — v u /A and I v = /4. The values E nu<nv denote the 

energies which are determined by the EBK-quantization conditions ( |33[ ) for either D- or iV-boundary 
conditions. For simplicity of notation we do not write an additional index for the boundary conditions. 
The quantities u u and v v are the Maslov indices. They are approximated by the uniform approximation 
given in ( p3| ) . After a change of variables from I u and I v to a and E one obtains 

OO .^2 

d(E)= V / do J(Q, J B)e 2 ™ (/ "-^ /4)+2 " m(7 "- iy " /4) , (38) 



n,m=— oo 



where J(a,E) is the Jacobian of the transformation. For values of a and E where a = a\E/(2c) is 
not small one can neglect the a-dependence of the Maslov phases and the Jacobian is given by 



J(a,E) 



dl u dl v dl v dl u 



1 

2E 



la 



Iv 



dlu 
da 



(39) 



dE da dE da 
which is a function of a only. 

We discuss in the following the different contributions to the integrals in (|38|). In doing this we 
assume that the eccentricity of the ellipse is not very small. If it is very small then the unstable and 
stable orbits along the two axis of the ellipse cannot be treated as isolated orbits, since they result 
from the break-up of a torus of the circular billiard. An investigation of the semiclassical contributions 



of these orbits for small eccentricity was carried out in [19| 



2.3.3 The area contribution 

The term with n = and m = in ( p8| ) and (^) is the only one with a non-oscillatory integrand. It 
has been shown in general that this term corresponds to the area term of Weyl's law @ 

/•OO /-OO A 

d , (E)= dl u dI v S(E-E nu , nv ) = — , (40) 

JO Jo 47T 

where A is the area of the billiard system. For a half-ellipse A = irab/2. 
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Figure 1: The curve (I v (a), I u (a)) for a half-ellipse with axis-ratio r = tt- The action I v is negative 
if a < — c 2 . The line from the origin marks the point corresponding to a = 0. 

2.3.4 Contributions of orbits with a > 



For all other terms the integrand in (^8|) is oscillatory and the main contribution comes from values of 
a where the phase of the exponential function is stationary. We consider first contributions from a > 
and approximate them by a stationary phase approximation. We neglect a possible o-dependence of 
the Maslov phase, i.e. we assume that the energy is large enough so that a = ayE/(2c) is not small 
at a stationary point and the Maslov index can be considered constant. Then the Jacobian is given by 
( |39|) and we drop its second argument since it is energy- independent. The integrals in ( |38[) can then 
be considered to be integrals over the energy surface E = const, of the billiard system. This energy 
surface is plotted in figure 1. The stationary points are those points in this figure where the slope of 
the curve is rational. The curvature of the energy surface is positive for a > and negative for a < 0. 
The stationary phase condition is 

d 

— {nl u + ml v ) = 0, (41) 

and it coincides with the condition (^) for periodic orbits. It has a solution for n = 3,4,... and 
1 < m < n/2 and for the corresponding negative values of m and n. The term for negative values 



(— n, — m) is the complex conjugate of the term for (n,m) as can be seen from Eq. (38) since J is 
real. For that reason we calculate the contributions for positive values of n and m and add at the 
end a complex conjugate contribution. The result of the stationary phase approximation for a term 
corresponding to a pair (n, m) is 

T( e , exp{27rm(J u - u -f) + 2mm(I v -f) + i\] 

J \ a n,m) r~aT T 7ttT~ ' ^ ' 

where a e n m is the value of a at the stationary point. The superscript e denotes that we consider in 
this section contributions of periodic orbits with elliptical caustic. 
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The second derivatives of the actions are given by 



d 2 L, k 2 / VEab 



da 2 4ac 2 \ w ir\fW - n 

d 2 i v k 2 



I v ■ (43) 



da 2 4ac 2 v 

This relation is also valid for a < 0. At the stationary points one further has 



] e FIT 
[ n ' m) ~ Ann^/E da ' 



which follows from Eqs. (|l~2|) and (|Ig|). 

Combining all results one obtains the joint contribution of the pairs (n,m) and (— n, —m) 



d e n ,m( E ) = V — i f== — cos ( kl n. m - T; n < - T; m < + t , (45) 

s irk I 2nab _ ]e "2 2 4 



where k = y/E is the wave number. The results for the two half-ellipses are identical since the Maslov 
indices are = 3 and v% = 4 for D, and u^ = 3 and v% = for iV and thus the argument of the cosine 
differs by a multiple of 2ir. 

For the derivation of Eq. ( [45| ) we assumed that the parameter a is large at a stationary point. If this 
is not the case then the approximation is more complicated. Then one has to take into account the in- 
dependence of the uniform approximation for the Maslov-phase in the exponent when one determines 
the stationary points. This leads to a replacement of condition ([II]) by 

d 

0^( nn u + mn v ) = , (46) 

where n u and n v are expressed as functions of a by Eqs. (p3|). The same modification applies also 
to tori with small negative values of a. Condition ( |46| ) does not correspond to the condition for 
periodic orbits ( |4l| ) any more, but in the limit E — > oo the previous stationary phase condition (|4l| ) 
is recovered. This is similar to the treatment of whispering gallery orbits in the circular billiard . 
(For contributions of whispering gallery orbits in the stadium billiard see |3^].) Their semiclassical 
description involves effective lengths which are different from the lengths of periodic orbits. We do 
not further discuss the modifications due to Eq. (f46|), a detailed investigation of contributions of tori 



near the separatrix is carried out in [37|. Finally, we remark that we did not consider corrections 
for whispering gallery orbits with a ~ b 2 in this section. These corrections go beyond the EBK- 
quantization 



2.3.5 Contributions of orbits with a < 

In the case of negative a the main contributions come again from the stationary points of the integrand 
in (|38|), and as before we assume the a is not small at these points. However, now the stationary points 
can also lie outside the integration range (when the corresponding torus is complex) or very close to the 
lower end of the integration range. For these cases a stationary phase approximation is not appropriate 
and one has to use a uniform approximation instead. We apply here the uniform approximation of 
Berry and Tabor § for integrals with a stationary point that can lie near the boundary or outside 
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the integration range. This approximation can be written in the form p2| 



da g{a) e 



ikf(a) _ g(a*)V2iri(3 ik f( a *) 
, j g(«o) -ifc/Cap) 



0(a* — ckq) + 



-sign(a* 



OO 1 



A' 2 



(47) 



where (3 = sign(/"(a*)), A = v /2A;/3(/(q ) — /(a*)), a* is determined by /'(a*) = 0, and 6(a) denotes 
the Heaviside theta function. The terms on the right-hand side of Eq. ( fl7| ) can be interpreted in the 
following way: the term multiplying the O-function is the stationary phase approximation of the 
integral, the last term is the contribution from the boundary (which can be obtained by an integration 
by parts), and the remaining term is an interference term between both. As will be seen in the following 
section the sum over all boundary contributions gives the Gutzwiller expression for the semiclassical 
contribution of the stable orbit along the minor axis of the billiard and its repetitions. When the ratio 
b/a of the two half-axis is decreased then new tori of periodic orbits arise through bifurcations of the 
stable orbit (and its repetitions). At the bifurcation point the new tori coincide with the stable orbit 
and have zero extension, and as b/a is decreased further they separate from the stable orbit. This 
is expressed by formula (|47|). Far away from the bifurcation the interference term can be neglected 
and the torus is well separated from the stable orbit. However, near the bifurcation, the contributions 
of the torus and the stable orbit cannot be separated, they rather give a joint contribution. At the 
bifurcation the Gutzwiller contribution of one repetition of the stable orbit diverges and this divergence 
is cancelled by a similar divergence of the interference term. 

We consider in this section only the stationary phase term and the interference term in Eq. (E7]). 
The sum over the boundary contributions is performed in the next section. The stationary phase 
condition is again given by Eq. (41) which has a solution for n = 3,4, .. . and 1 < m < n/2 and for 
the corresponding negative values of m and n (since we consider a half-ellipse n does not have to be 
even). As before we calculate only the contributions for positive values of n and m and add at the 
end a complex conjugate contribution. 

The second derivatives of the actions are given by Eq. (^) and the value of the Jacobian at a 
stationary point by Eq. (^J|) with the difference that now the quantities are given a superscript h to 
signify the contributions of tori with a hyperbolical caustic. The value of (3 in ( [47| ) is —1 since the 
curvature of the energy surface is negative for negative values of a. Inserting these values into (PTJ 
one obtains 



'2o2 



] h F ( 



2 > 



771 

x 'n . m 



) 



irk 



mr 



2nab 



n.m 



h )2 



0(aS, m + c 2 ) cos(fcZ* - ^-nv 



+ 



si g n «,m + c 

2^ 



dX 



2 

sin(A;Z 



^ h 



1^ 
4' 



h 

n,m 



X 2 /2 



X 2 



(48) 



where A = ^2k(l\, 
= 2 for N. 



2nb) and the Maslov indices are v£ = 4 and v% = 2 for D, and 1/^ = 2 and 



2.3.6 The contribution of the stable orbit 



We sum now over all contributions from the boundary a = —c of the integrals in fl38|). These 
boundary contributions are given by the last term in (^) and they exist for all values of n and m, 
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even when there is no stationary point for these values. The only exception is when both values are 
equal to zero because then the integrand is non-oscillatory and yields no boundary contribution. It 
has been shown in general by Richens that the contribution of a stable periodic orbit can be obtained 
by summing over all boundary contributions in an integrable system |22|] . 
The values of the actions and their derivatives at a = — c 2 are given by 



by/E dl u 



Iu\a=—c 2 ., 

7r oa 
II 2 -0 ^ 

1 v\a=-c I — U n 

oa 



E . b 
arcsm — 



2ttc a 
E 



(49) 



. - c 2 2c 

and J{— c 2 ) = 6/(47rc). Summing over all boundary contributions one obtains 

y, , ib exp{2inkb - *fnv* - l fmv%} 
n m= _ 0O 47rfc irm-n arcsin £ 

g b sin(2nfc6-fn^) (5q) 
-, irk 2 sin(n arcsin -) 

n=l v a y 

where the prime at the first sum indicates that the term (n, m) = (0, 0) is omitted. Furthermore, the 
relation 

- HT 1 

L — ' z — Tim sin z 

m=—oo y ' 

has been used [35| and the fact that v% = 2 for D and N. 

The semiclassical contribution of the stable orbit and its repetitions d s {E) is divergent when 
n arcsin (6/a) = mi:. This coincides with the condition Eq. ( p2[ ) for the appearance of a new torus. 

2.3.7 The perimeter contribution 

At the boundary of the integration a = b 2 one obtains also a contribution to the energy density. It is 
obtained by the approximation 

boundary contribution of ( H da g(a) e ikf{a) \ » _lii"lL*fc/(aa) , (52) 

U-oo J k J {ai) 

which follows from an integration by parts. Near a = b 2 the actions are expanded in powers of 
e = b 2 - a: 



In 



3nab 



2a\/E? ,7r c, ^/E^.tt c, . , 

4 « E -,- -e F (— , -) . 53 

7r 2 a na 2 a 

It follows that J(a) ~ —I v (2E)~ 1 dI u /da and it's leading term is proportional to e 1 ! 2 for small e. For 
terms with m ^ the exponent in the integral (^) depends in leading order linearly on e and the 
corresponding boundary contribution vanishes. The contributions of the terms with m = and 
are given by 

00 rb 2 / J v QJ u \ 27T 
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exp{2mnl u - —nv u ) 



n=l 
oo 



2lTnE 



=6 2 



2ttE 



E 



sin(fn<) 



n=l 

aE (c/a) 
47rA: 



(54) 



where = 3 has been inserted for Z) and N. The term ( |54"| ) is the perimeter contribution —L/(8irk) 
of the outer arc of the half-ellipse to the mean spectral density, where L = 2aE(e). 



2.3.8 The contribution of the unstable orbit 

We now consider the question how the contribution of the unstable periodic orbit can be obtained. 
This orbit has an a-value of zero, and at this point the original EBK-quantization is discontinuous 
and inaccurate. It has been first pointed out by Bogomolny that the semiclassical contribution of the 
unstable orbit can be obtained by using the uniform approximation for the semiclassical phase near 
the separatrix [40|. We follow his method in this section. 

In the context of the uniform approximation it is more convenient to work with the a-variable 
instead of the a-variable. For that reason we start again with equation ( |37| ) and change the integration 
variable from I u and /„ to a and E and obtain 



d N ' D (E)= / d « J («' E ) exp{27rin(/ u 



PsA a )< 



TT 



+ 2irim(I v 



TT 



)} , (55) 



where the new Jacobian is given by 

dn u dn v dn v dn u 



J(a,E) 



dE da dE da 
din /dl v | 2 dPsA~ 
dE \da it da 



(56) 



a, 



dh fdlu 
dE V da 



7r da 



We evaluate now the contributions to the integrals in fl55[ ) from the vicinity of a = 0, and we restrict 
to contributions with n = 2m since these terms are the ones which correspond to the unstable orbit 



as has been discussed in section 2.1 



The derivative of PsA a ) with respect to a is given by 



^=PsA*) 



7T 



1 ,,1,1 

. ■ H — log \a\ ty( — h ia) 

4cosh(vra) 2 1 1 4 y 2 ' 



1 A 

4 v 2 



ia 



(57) 



where *&(z) is the logarithmic derivative of the T-function [35]. For small values of a one further has 

^E{a - c) 

~ 2^~ 

Zc\fE 



IT 

2VEc 



a a 8c\^E a a + c 

7, 7T lo £ i-i + - log ' 

ZTT ZTT \a\ 



b 



TT 



a a , 
+ - + - log ■ 

7T 7T 



a 



(58) 



With these equations the Jacobian J is evaluated near a = 
1 fa. 2c . a + c . a 



J(a,E) 



2nVE Vvr 



a 

- log(8c 



2c , a + c 
E --log-— ±- 

tt o 2cosn(7raJ 



a TT /I --\ T A 

-M- +l a ) + n- 2 



ia 



■ (59) 
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In expression (58) for J we have set all appearing a-values equal to zero except in the argument of 
the cosh- and ^-function. The reason for this will become clear in the following. Inserting Eqs. fl58|) 
and fl59|) into (]55|) one obtains 



d»> D (E) = 2Re V f da t^l ( 1 i og (8 C v^) - * log i±£ ± fL 
u v ' £± X J 2 r K\[E \7T v 7 vr & 6 2cosh(vra) 

_ 7T-[*(^ + is ) + - ex P \^imaVE + 4imalog \ , (60) 
27r 2 2 / t b J 

where the relations fls{a) + Ps{—<x) = 7r /4 and /3a(<5) + Pa{—ol) = 37r/4 have been used in the 
exponent. 

As one can see, the exponent in (^) depends linearly on a and for that reason one would not 
expect to obtain a contribution of the order of an unstable orbit from a = 0. The reason that one 
does obtain this contribution nevertheless is that the Jacobian in (|59|) has poles with Re a = 0. More 
precisely, ^(1/2 + ia) has poles at i(l + 1/2) with residuum i for I = 0, 1, . . ., ^(1/2 — ia) has poles at 
— i(l + 1/2) with residuum —i for Z = 0, 1, ... , and 1/ cosh(7ra) has poles at i(l + 1/2) with residuum 
— £(— l)*/7r for I £ The contribution of the unstable orbit is obtained by closing the integration 
contour in the upper half plane and using Cauchy's theorem 

OO ( OO / -i \yxi - - 

d^ D (E) = 2Re V I Y a[ ~ exp{4imak - 2m log - 4m/ log ^} 



CO 



a(-l) m (-l)' a + c a + c.} 

± > expj 4imak — 2m log — 4mi log \ > 

^ 2vr/c Pl B 6 B b J j 

J a ( — l) m exp{4imafc} a(— l) m exp{4imafc} 1 

2 e ^ 1 j 27Tfc (2+£)2m _ ( a+c)- 2m ± 2 ^ ( a+c )2m + ( a+c)-2 m j 
oo r 

= } — cos(4ma/c — irm) < — — — ± — — — > . (61) 

_^ irk 12 sinh(2m arcsinh|) 2 cosh(2m arcsmh|) J 

This is the contribution of the unstable orbit which runs along the x-axis that is part of the boundary 
of the half-ellipse. 

2.3.9 The joint contribution 

Combining all results of the previous sections and adding the contributions for D and iV one obtains 
the trace formula for the full ellipse. We give it here for the level density d(k) = 2kd{E) in terms of 
the wave number k 



rn-li 



d(k) = d(k) + )^ 2^ V ' =cos(kl - —n + - 



E E 



/ r~i 7 7 Alh V (— — 



©("n,m + c 2 ) cos(kl^ m - vrm - ^) 



sign(<, m + c 2 ) sin(M* m - X 2 /2 - vrm) 



2¥ 7a 



^ 45 sin(4m/c&) ^ 4a cos(4m/ca — m7r) 

^ vr 2sin(2m arcsin^) ^ vr 2 sinh(2m arcsinhf ) ' 
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where d(k) is the mean spectral density and A = y 2k(l\ m — 2nb). 

The same method can be used to obtain also the trace formula for quarter ellipses with the four 
different boundary conditions given in section 2f2. We give here only the result. The modified EBK- 
quantization conditions are 



DD : 
DN 
ND 
NN 



In 

Iu 
Iu 

I, 



j_ 1 , /3 A (a) I v 1 A (-a) 
n u + - H — = n v + - H 

I Z 7T 



2 7T 
^ 2 7T 

2 TV 



y = ™« + — - — 

4 ,1 Ps(-a) 

-2 =n " + 2 + —r- 
h 

y =n,+ 



7T 



'ifl! = 1,2, 
T^Ui Tlv — 1)2, 

; — 1)2, 



(63) 



where I v /2 with definition (10) is the value of the action in a quarter ellipse. The trace formula for a 
quarter ellipse is given by 



d(k) 



n-l 

OO 2 



2ka e j %m 



je -p. /7T e \ 



n—2 m=^ 



7T 



n7r 



2nab 



u n,m 



: COs(/c/^ , 



7T TT^ 

—nu,, — irmv,, + — , 



n-l 
OO 2 

+ E E 

n=2 m= l 



2fc|a^ r 



/ft y(E } \ 



7T 



nvr /( 



2nab 



^/b 2 -a h 



Jh \( K h 
n.mJ \ n; 



+ 



sign(a^, m + c 2 ) r<x> sin(fc# _ - X 2 /2 



dX 



X 2 



b I sin(2m/c6 - § mv%) , cos(2mkb - § mz^ - f - 1)) 



m= 
oo 



'i 7T 2sin(m arcsin-) 



+ 



+ E 

m=l 



2 cos(m arcsin -) 



cos(2mka — ^m(y u + u v )) cos(2mka — ^m(u u + v v ) — \(v^ — 2)) 



2 sinh(m arcsinh 



2 cosh(m arcsinh 



(64) 



where m runs now over half-integers in the sums over tori, A = J2k{l^ m — 2nb), and the Maslov 
indices are given by 

k 



DD : 


K 


= 3 


< 


= 4 


DN : 




= 3 


< 


= 2 


ND : 


< 


= 3 




= 2 


NN : 




= 3 


< 


= 



, h — 4 

K = 4 

4 = 2 

4 = 2 



^ = 3 
<4 



1 

3 

1 . 



(65) 



The sum (v u + z/y) in the last line of Eq. ( pi}) can be taken either at a > or a < since it is invariant. 

We remark that in Eqs. (^) and (p4|) we did not include modifications for tori close to the separatrix 
that were discussed in section |2.3-4 



2.4 Numerical results 



In this section we compare the semiclassical approximation of the last section with exact quantum 
results. For this purpose we consider a Fourier cosine transform of the oscillatory part of the level 
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Figure 2: The Fourier cosine transforms D QM (x) (full) and D sc (dashed) of the level density, calculated 
from quantum energies and periodic orbits, respectively. 
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Figure 3: The differences D QM (x) — D sc (x) and D EBK (x) — D sc (x) between Fourier cosine transforms 
of the level density. 
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Figure 4: The differences D QM {x) — D sc (x) between Fourier cosine transforms of the level density 
with (full) and without (dashed) contributions from interference terms to D sc (x). 



density with a Gaussian damping term and a cut-off 



D(x) 



1 



v max JO 



dk [d(k) — d(k)] cos(fcr) exp(— z 



U2 

""max 



(66) 



This function has peaks in the vicinity of the lengths of periodic orbits. 

The calculations were carried out for a quarter ellipse with a = 1.1 and b = 0.9 and with Dirichlet 
boundary conditions on all sides. The damping factor was chosen to be z = 5. The function D(x) 
was evaluated in three different ways and the results are denoted D QM (x), D sc (x) and D EBK (x). The 
first function D QM (x) was obtained from the quantum energies. These energy levels were determined 
up to energy E max = 20 000 by solving the two Mathieu equations (E7|) and fl2~g|) numerically, which 



corresponds to 1 198 energy levels. The second function D sc (x) was determined by using Eq. (64) with 
all periodic orbits up to length I = 12 and n = 30, and the third function D EBK (x) was determined 
from the semiclassical energy levels that are solutions of the EBK-conditions (|63|). 

In figure 2 the results for the quantum spectrum and the trace formula, D QM (x) and D sc (x), are 
compared. Both curves are in good agreement and can hardly be distinguished. For that reason we 
plot the difference between both curves in figure 3 (dashed line). One can see that the semiclassical 
error is approximately one order of magnitude smaller than the function D(x). Figure 3 shows also 
the difference between D EBK {x) and D sc {x) (full line) which is much smaller than the difference 
between D QM (x) and D sc {x). This shows that the error which was introduced by deriving the trace 
formula from the EBK-energies with stationary phase and uniform approximations is much smaller 
than the original error of the EBK-quantization. This cannot necessarily be expected because the EBK- 
quantization condition, the trace formula and the stationary phase approximations which connect both 
approximations are all only valid in leading order of ft. A similar result was observed previously for 
a circular billiard with a singular magnetic flux line [41]. The smallness of D EBK (x) — D sc (x) shows 
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also that the modifications for contributions of tori near the separatrix due to Eq. ([!(]) are not large 
in the range where the numerical examination was carried out. 

We also show that it is important to include complex orbits, i.e. the interference terms, in (|G4|). 
In figure 4 the difference D QM {x) — D sc {x) is plotted where D sc (x) was calculated once with (full 
line) and once without (dashed line) the interference terms. One can see that the semiclassical error 
is much bigger in the second case, especially near x = 9. This is due to a torus of complex orbits with 
n = 5 and m = 1.5 which is close to becoming real. It becomes real for b/a ~ 0.809 which is near the 
present value of b/a = 11/9 « 0.818. 



3 The oval billiard 

We discuss now how the semiclassical approximation for the elliptical billiard has to be modified when 
the ellipse experiences a small perturbation which makes the system non-integrable. We consider in 
particular a perturbation which consists in a deformation of the ellipse into an oval-shaped billiard 



system introduced by Berry [3C]. It has a parameterization which expresses the radius of curvature R 



of the boundary as a function of the angle Hf between the tangent vector and the x-direction 

R{^) = l + <5cos(2$) , (67) 
where ^ lies in the range [0, 2n). From this definition the dependence of the x- and y-coordinates on 



^ follows as [30] 



s s 

xtfy) = sin(^) + - sin(^) + - sin(3^) 
2 6 

5 5 

y(^) = -cos(^) + -cos(^) - -cos(3tf) . (68) 
2 6 

The oval has an area A = ir — ir5 2 /6, a perimeter L = 2tt, two half-axes of lengths a = 1 + 5/3 and 
b = 1 — 5/3, and it is a deformation of an ellipse of order 5 2 . The billiard can as well be considered as 
a perturbation of the circle. Then the deformation is of order 6. 

Numerical examinations were performed for the deformation parameter 5 = 0.3 for which the half- 
axes of the oval billiard have the same lengths as those of the elliptical billiard in the last section. We 
consider again a desymmetrized version of the billiard which consists of a quarter oval with Dirichlet 
boundary conditions on all sides. For this system the shortest periodic orbits were determined. Figure 5 
represents all periodic orbits which are reflected up to five times on the oval-shaped boundary of the 
quarter oval. We choose to present them by corresponding orbits in the full oval since their structure 
can be seen more clearly there. All these orbits appear in pairs which are labelled by letters. With 
the exception of the pair a all pairs of orbits result from the break-up of tori of the ellipse. Since the 
perturbation is small, no subsequent bifurcations of orbits occur in this length regime. 

The differences in the lengths of orbits within a pair are shown in Table 1. With the exception 
of the first pair the length differences are very small. For that reason, the orbits cannot be treated 
as isolated orbits in a semiclassical approximation if the energy is not rather high. We discuss in 
the following the different uniform approximations by which the contributions of the pairs to the 
semiclassical trace formula can be treated. Most of the pairs (b, c, d, /, g, i,j, k, I) result from a break- 
up of tori of the elliptical billiard with a > 0, i. e. from tori which have a confocal ellipse as caustic. 
Their contribution can be treated by the uniform approximation for the generic break-up of tori, that 
will be described in the next section. 

Other pairs (e, h, m) result from the break-up of tori of the ellipse with a < 0, i. e. from tori with a 
confocal hyperbola as caustic. These tori all arose from a bifurcation of the stable orbit in the ellipse, 
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Figure 5: Periodic orbits in the oval billiard for S = 0.3. The desymmetrized versions of these orbits in 
the quarter oval represent all orbits which are reflected up to five times on the oval-shaped boundary. 
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At 
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3.00259 
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-0.00035 


h 
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7.49347 
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0.02427 
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7.54064 
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-0.01176 
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3.09101 
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0.00001 
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5.88609 


0.00007 


0.00031 


e 


k 


8.12883 


0.00029 


0.00288 


e 
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9.70973 


0.00149 


0.07798 


e 


m 


y. I 2L6Z 


0.00135 


— U.UD491 
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2b 


5.67096 


0.01411 


0.61722 


e 


3b 


8.50644 


0.02116 


3.17563 


e 



Table 1: Properties of pairs of orbits in the quarter oval. The different columns show the mean length I 
the length difference Al and the quantity At = Tr Mi +Tr Mi — 4 for a pair of orbits. The type-column 
specifies whether the pair arose from the break-up of a torus of the ellipse with a confocal ellipse or a 
confocal hyperbola as caustic, respectively. 

as the eccentricity of the ellipse was increased starting from the circle. These tori are not necessarily 
well separated from the stable orbit in the ellipse, as is expressed by the interference term in (|4S|), and 
thus the usual formula for the break-up of tori cannot be used. 

The correct contributions of these pairs to the trace formula can be obtained by considering the 
oval not as a deformation of the ellipse but as a deformation of the circle. Then the pairs (e, h, m) 
are obtained not from the break-up of tori but from bifurcations of the stable orbit along the vertical 
axis of the billiard (the first orbit of the pair a in figure 5). This stable orbit will be denoted by 
a s in the following. In particular, the pair e arises from a bifurcation of the 3- fold repetition of 
the orbit a s at 5bif = 3/13 ~ 0.2308, the pair h from a bifurcation of the 4-fold repetition of a s at 
5 Hf = (39 - 24\/2)/41 ra 0.1234, and the pair m from a bifurcation of the 5-fold repetition of a s at 
Sfoj = (63 — 24\/5)/121 ~ 0.0771. For the pair e the bifurcation is shown in figure 6 by plotting the 
orbits for values of S near 5uf ~ 0.2308. This bifurcation results in four new orbits in the full oval, 
the right orbit in figure 6 which can be traversed in both directions and thus represents two orbits, 
and the left orbit and its mirror image (obtained by reflection on the y-axis) which is not plotted. 

It follows from the above considerations that the contributions of the pairs e, h and m are described 
by formulas for bifurcations that are discussed in section [3.2| . This is a very general property. If one 
considers an integrable system which depends on an external parameter and in which bifurcations occur 
that result in the appearance of new tori, then in a slightly perturbed system there are bifurcations 
which result in a finite number of new orbits. Thus the break-up of tori which are close to a bifurcation 
is described by formulas for bifurcations in which a finite number of orbits is involved. 
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Figure 6: The period-tripling bifurcation which yields the orbits e of figure 5. The orbits are shown for 
5 = 0.230, 5 = 0.235 and 5 = 0.240 (from below), and the bifurcation occurs at 5uf = 3/13 0.231. 



3.1 Broken torus contribution 

If an integrable system is perturbed then all tori of periodic orbits break up and only a finite number 
of periodic orbits remain (except for degenerate cases). According to the Poincare-Birkhoff theorem 
the number of remaining orbits is even, half them are stable and half of them are unstable. The 
Poincare-Birkhoff theorem does not specify the total number of remaining orbits, but generally this 
number is two [42]. For the semiclassical treatment of this case, i. e. the break-up of a torus into one 
stable and one unstable orbit, a uniform approximation was derived by Tomsovic, Grinberg and Ullmo 
This approximation extends previous results of Ozorio de Almeida [18, 43]. The formula can 
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also be applied to cases in which, due to the presence of discrete symmetries, the break-up of a torus 
results in n stable and n unstable orbits with respectively identical actions, periods, stabilities and 
Maslov indices. The uniform approximation for the contributions of these orbits to the level density 
can be written in the following form 
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(69) 



In case of the break-up of a torus into 2n orbits due to the presence of discrete symmetries the 
equation has to be multiplied by n. The indices u and s correspond to the unstable and stable orbit, 
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respectively. Furthermore 

s= Su + s 1 AS 



T u +T s 



AT 



(70) 



and S, T, u, M and r denote the action, period, Maslov index, stability matrix and repetition number 
of an orbit, respectively. 

In the following we will use a slightly simplified version of (^9|) which is obtained by applying the 
relation ^2(2) = 2 3±(z)/z — Jq(z). The term 2 Ji(z)/z yields a modification of the prefactor of the 
Ji-Bessel function which is one order of h smaller than the previous prefactor of Ji in (|69|). Since we 
consider only the leading order semiclassical approximation we neglect this term and obtain 
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(72) 



This approximation, as well as (|69|), interpolates over the whole regime from the torus contribution 
to the contributions of isolated periodic orbits in the Gutzwiller approximation. In the limit AS — > 
the mean amplitude A diverges, but the product Ay/\ AS\ remains finite, and the whole expression 
yields the Berry- Tabor term for the semiclassical contribution of a torus. In the opposite limit where 
AS/% >> 1 one can replace the Bessel functions by their leading order term for large arguments and 
obtains the Gutzwiller approximation for the contributions of isolated orbits. Finally we note that 
equation ( |7"l| ) in combination with the definitions ( |7"0| ) and ( [72] ) is invariant under exchange of the 
labels u and s. Thus the definitions in (f7C|) and (72) can also be formulated in terms of an orbit 1 and 
an orbit 2 without specifying which of them is stable and which unstable. 



3.2 Contribution of bifurcating orbits 

There is only a finite number of generic forms in which periodic orbits in two-dimensional systems 
bifurcate. They are characterized by normal forms that describe the classical motion in the vicinity of 
a bifurcation. The different normal forms for generic bifurcations of periodic orbits in two-dimensional 
systems were derived by Meyer [p]] and Bruno [12, 13] and are discussed in [43, H. Altogether there 
are six different kinds of bifurcations which are classified according to the lowest repetition number 
m of a central orbit for which the bifurcation occurs. (A special case is m = 1 where there is no orbit 
before the bifurcation.) There is one kind of bifurcation for m = 1 up to m = 3, respectively, two kinds 
for m = 4, and one for m > 4. The corresponding bifurcations are period-m-tupling bifurcations, i. e. 
the primitive periods of the arising periodic orbits are m times the primitive period of the bifurcating 
central periodic orbit at the bifurcation. 

We restrict the discussion in the following to the case m > 4. This kind of bifurcation is the first 
that can be observed when an integrable system is perturbed since the higher repetitions of a stable 
orbit bifurcate earlier than the lower repetitions. It has the following form: a central stable periodic 
orbit bifurcates and two new periodic orbits appear, one stable and one unstable, while the central 
orbit remains stable. The new appearing periodic orbits are called satellite orbits. 



For the case m > 4 a uniform approximation was derived in [17| which interpolates over the regime 
from the bifurcation up to regions where the orbits can be described in the Gutzwiller approximation. 
This uniform approximation was obtained by an extension of the method of Ozorio de Almeida and 
Hannay [14|. 
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The final formula in fll7| was given for bifurcating periodic orbits in billiard systems and for 
bifurcating periodic orbits without turning points in systems with potentials. The derivation in jl?]] 
can also be applied to orbits with turning points, since only the Maslov index has to be changed. 
The following uniform formula for the contributions of bifurcating periodic orbits to the level density 
includes all these cases 



, 9(e) 2tt\AS\ Jt /AS\ S tt \ 

*^ = -^v n { AJo UJ cos U-rJ +AAJi vn 

. A [\AS\ r A v cos (I + |X 2 -|(^-/3)' 
+sign(£) — Wi—i / dX 



tt/j V h J A X 2 -Kfimr y/\Tr M 




Here A = y 2/3[So — S]/h and £ is a parameter that is positive when all three orbits are real, and 
negative when only the central orbit is real. The values of u and (5 can be obtained from the Maslov 
indices of the (real) periodic orbits, since u$ = u + sign(e) • /3, v u = v and u s = u — j3. The other 
quantities are defined in ( |70| ) and (f72|). The repetition number of the central stable orbit is mr where 
m is the lowest repetition number for which the bifurcation occurs. (When the m-th repetition of 
a periodic orbit bifurcates, then the bifurcation occurs also for all repetition numbers which are a 



multiple of m.) Eq. (73) is valid as long as no subsequent bifurcations of the participating periodic 
orbits occur before they can be considered isolated. Note that formula (^) in combination with 
definitions @) and © is again invariant under exchange of the indices u and s. The different terms 
in Eq. (|7^) can be interpreted in the following way. The first term with the two Bessel functions is 
a joint contribution of the two satellite orbits. It has a form which is identical to the broken torus 
contribution J7l|). The last term in (ff^) is the semiclassical contribution of the central stable orbit, 
and the remaining term is an interference term between satellite orbits and central orbit. 

The properties of the three orbits that participate in the bifurcation are summarized in Table 2. 



= 1 


£ < : £o stable, uq = v — 1 

£ > : £o stable, uq = v + 1 
£ u unstable, v u = v 
£ s stable, u s = v — 1 
Sq > S u > S s 


p = -l 


£ < : £o stable, uq = u + 1 

£ > : £o stable, uq = u — 1 
£ u unstable, u u = u 
£ s stable, u s = u + 1 
S s > S u > Sq 



Table 2: Properties of orbits which participate in a generic period m-tupling bifurcation with m > 5. 
£o denotes the central orbit, and £ M and £ s denote the unstable and stable satellite orbits, respectively. 

The bifurcations in the oval billiard that lead to the pairs e, h and m are not generic. This is 
due to the symmetries of the billiard system. But the normal forms which describe these bifurcations 
agree with normal forms of generic bifurcations, namely with those for double the repetition number. 
For example, the pair e results from a bifurcation of the 3- fold repetition of the orbit a s , but its 
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normal form corresponds to that of a generic bifurcation with repetition number m = 6. This follows 



from the treatment of bifurcations in the presence of symmetries [45, 46]. It can be understood by 
considering a Poincare section of surface perpendicular to the central orbit. For a generic period-m- 
tupling bifurcation with m > 4, the map from the Poincare section of surface at some starting point 
to the Poincare section of surface after m traversals of the central orbit has 2m fixed points near 
the central orbit (after the bifurcation); m of them correspond to the new stable orbit and m to the 
new unstable orbit, since both orbits cross the Poincare section of surface m times before they close. 
For the considered bifurcations in the oval billiard there are 4m instead of 2m fixed points in the 
vicinity of the central orbit, since two new stable and two new unstable orbits arise, but the number 
and arrangement of stable and unstable fixed points is the same as for a generic period-(2m)-tupling 
bifurcation. 

As a consequence, also the period-tripling and period-quadrupling bifurcations of the stable vertical 
orbit a s in the oval billiard can be described by the formula for generic bifurcations with m > 4. The 
only difference is, that m has to be replaced by 2m in Eq. (f7^) and the whole formula has to be 
multiplied by two, if one considers the bifurcation in the full oval. 

In the desymmetrized billiard, the quarter oval, the orbit a s runs along a part of the boundary of 
the billiard and Eq. ( |73| ) has to be modified in a different way. We describe this by first looking at 
the Gutzwiller contribution of the n-th repetition of the orbit a s which is modified to take account of 
the fact that a s is a boundary orbit 
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(74) 



where b = 1 — 5/3 and v = arcsin(i/(l — 5/3) /(l + 6)). The second line in Eq. (f73) is given in 
dimensionless units. The quantities with index 1 are those of the vertical stable orbit in the half oval 
with Dirichlet boundary conditions on the x-axis. Furthermore, T 2 = Ti, S 2 = S\, M 2 = —Mi and 
^2 = v\ + 1 + 2([nu/7r]mod2). Eq. (|74| ) has the form of a sum of contributions from two orbits, each 
with half the usual amplitude. Bifurcations of the n-th repetition of the orbit a s occur when 5^0 
and one of the two terms in (]74] ) diverges. The approximation ( [74]) is only valid if the orbit is well 
separated from a bifurcation. Near a bifurcation it has to be modified by replacing the term which 
diverges at the bifurcation by Eq. ([7^) where m is twice the lowest repetition number for which the 
bifurcation occurs and r = 2n/m. If both terms are close to a bifurcation then both terms have to be 
replaced. 

In the following we apply the uniform approximation (f/"3|) and its modifications in order to describe 
the break-up of tori of the ellipse with a < 0. A comparison with the semiclassical contributions of the 
tori in the elliptical billiard (section |2.3.5| ) shows that this effectively corresponds to an application 
of the broken torus approximation ( |7i~D to the torus term in Eq. (f4"8|). The interference term and the 
contribution of the bifurcating stable orbit do not change their form. This is related to the fact, that for 
a generic period-m-tupling bifurcation with m > 4 the differences of the actions of the satellite orbits 
grow more slowly than the differences between the actions of the satellite orbits and the action of the 
central stable orbit. This is not the case for all bifurcations. For example, there is a period-doubling 
bifurcation in the ellipse for a ratio of the half-axis a/b = \[2. This bifurcation is again described by 
the approximation Eq. ( [48|) (plus the contribution of the stable orbit), but in this case it is not correct 
to apply the broken torus approximation ([71]) to the torus term in Eq. ( [48D in order to describe the 
break up of the torus. Instead one has to use the uniform approximation for a generic bifurcation 
with m = 4 which has a more complicated form and is described by a diffraction catastrophe integral 
for the catastrophe Xg |j"41 . 
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3.3 Numerical results 



In the following we examine different approximations for the spectral density numerically. For that 
purpose the energy levels of the quarter oval with 5 = 0.3 were determined by a boundary element 
method up to energy E max = 20 000. This corresponds to 1190 different energy levels. 

In order to compare semiclassical results with quantum mechanical results we consider again the 
Fourier cosine transform D{x) with a cut-off that is defined in (^). The damping factor is again z = 5. 
Figure 7 compares determinations of D{x) from the quantum mechanical spectrum D QM (x) and the 
uniform approximation D C ua(x). For the calculation of D CUA {x) all periodic orbits up to length I = 12 
and up to 30 reflections on the boundary were included. Both curves are in good agreement and it 
is difficult to see the differences. For that reason we plot in the following figures the semiclassical 
error, i.e. the difference between the quantum mechanical curve D QM (x) and different semiclassical 
(or uniform) approximations for it. 

The most basic semiclassical approximation for a quasi-integrable system, i. e. for a slightly per- 
turbed integrable system, is the quasi-torus approximation. In this approximation the pairs of periodic 
orbits are treated as if they still contribute like a torus of periodic orbits. We use here a quasi-torus 
approximation in which the semiclassical amplitudes and phases are determined directly from the 
corresponding periodic orbits. This is done by setting AS equal to zero in Eq. (|7"H). This approxima- 
tion is used for all pairs of orbits except for the pair a and its repetitions. For this pair the action 
difference is already quite large, and the orbits are treated as separate periodic orbits. In more detail, 
the contribution of the n-th repetition of the unstable orbit a u is approximated by 



cos(2nka — n^) cos(2nka — n 2 



2sinh(nu) 2cosh(nu) 



(75) 



with a = 1 + 5/3 and u = arccosl^^l + <5/3)/(l — 5)), and the contributions of the stable orbit a s 
and its repetitions by Eq. (f74[). 

The semiclassical error for the quasi-torus approximation D QM (x) — D QT (x) is plotted as dotted 
line in figure 8. This approximation already works relatively good. This is due to the fact that the 
action differences within pairs of orbits are still relatively small and that the orbits are not very close 
to a bifurcation. However, in comparison with the semiclassical error for the ellipse in figure 3 the 
semiclassical error for the quasi-torus approximation in figure 8 is clearly larger. The scales of the 
plots differ by approximately 50%. In order to improve the approximation we apply Eq. (^) for the 
pairs of orbits which arose from a bifurcation, i. e. for the pairs e, h and m, but we still set the action 
differences in formulas ( |7l| ) and ([F|) equal to zero. Note that with this substitution Eq. ( f73|) ha s 
exactly the same form as that of the contribution of a torus near a bifurcation (see section ^ .3. 5| ). 
This approximation can thus be considered as an improved quasi-torus approximation which takes 
into account cases in which quasi-tori are close to a bifurcation. The corresponding semiclassical 
error D QM (x) — D QTI (x) is plotted in figure 8 as full line. One can see an clear improvement near 
x ~ 5.4. This length corresponds to the pair e which is that pair out of e, h and m which is closest 
to a bifurcation. This can be seen from the differences between the mean lengths of the pairs and 
the m-th repetition of the central orbit from which they arose by a bifurcation: l e — 3l ag = 0.03636, 
lh — 4:l as = 0.34064, and l m — 5l aa = 0.72132. The pairs h and m are already well separated from the 
orbit a s and thus the inclusion of the interference term does not improve the approximation. 

For a further improvement of the semiclassical approximation we apply now the uniform approx- 
imations ([7l|) and (73) without setting AS equal to zero. The result D QM (x) — D UA (x) is shown in 



figure 9 as full line. For a comparison with the previous approximation the function D QM (x) — D QTI (x) 
is plotted again as dotted line. One can observe a clear improvement of the approximation at three 
places which correspond to the first three repetitions of the pair b. Table 1 shows that this is the pair 
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Figure 7: The Fourier cosine transforms D QM {x) (full) and D CUA (dashed) of the level density, calcu- 
lated from quantum energies and periodic orbits, respectively. 



< 




Figure 8: The semiclassical errors D QM (x) — D QT (x) (dotted) and D QM (x) — D QTI (x) (full) for the 
quasi-torus approximation with and without interference terms, respectively. 
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for which the length difference is largest, and the size of the semiclassical error increases with the size 
of Al. We note that for the decrease of the semiclassical error it is important to use the whole formula 
( |7l"| ) and not only the first term with the Jo-Bessel function. A criterion for the importance of the 
second term is the size of At in Table 1 which increases rapidly with the repetition number of the pair 
b. 

Even for the uniform approximation there is still a relatively large error at x ~ 9. This error is due 
to the fact that up to now we considered only contributions from real orbits to the spectral density. But 
before a bifurcation there are also contributions from complex orbits as is expressed by the interference 
term in Eq. (ff^) . The error at x ~ 9 is due to fact that the 5- fold repetition of the orbit a s is close 
to a second bifurcation which occurs at S^f = (75 — 24\/5)/61 ~ 0.3497. In a last step we improve 
the approximation by taking into account the contribution of these complex orbits. This is done by 
determining these orbits for values 5 > 5bif where they are real. From these data the properties of the 
orbits can be extrapolated into the complex region. The final result D QM {x) — D CUA {x) is shown in 
figure 10 as full line. For a comparison with the previous approximation the function D QM (x) — D UA (x) 
is plotted again as dotted line. As one can see, the error is clearly decreased near x = 9 by including 
the contribution of the complex orbits, and the final semiclassical error is of the same size as in 
the elliptical billiard. This shows that by using the uniform approximations ( |7l| ) and (|7^) the same 
accuracy of the semiclassical approximation can be achieved as in the unperturbed integrable system. 



4 Conclusions 

Most applications of semiclassical trace formulas in the literature have concentrated on cases where the 
periodic orbits are either isolated or occur in families. These approximation are valid only in restricted 
classes of systems. In general systems, it is common that semiclassical contributions of periodic orbits 
cannot be considered isolated from contributions of other periodic orbits in their neighbourhood. This 
requires the use of uniform approximations which take account of the underlying classical structures 
and yield joint contributions of neighbouring periodic orbits. In systems where the classical phase 
space has regular and chaotic regions this is a generic situation, and semiclassical quantization rules 
in terms of periodic orbits in these systems always involve uniform approximations. But uniform 
approximations can also be necessary in integrable or chaotic systems, for example when bifurcations 
of periodic orbits occur in integrable systems or when periodic orbits are close to creeping or diffractive 



orbits in chaotic billiard systems [47, 48]. 



In the first part of this paper we examined an integrable system, a billiard in form of an ellipse, 
and derived a trace formula for it. The elliptical billiard is a standard example for an integrable 
system, yet its semiclassical trace formula is more complicated than a summation over semiclassical 
contributions of isolated tori. These complications are due to the presence of the separatrix and of 
the stable periodic orbit. We were particularly interested in the treatment of the bifurcations of the 
stable orbit and its connection to generic bifurcations when the ellipse is perturbed. We derived a 
semiclassical trace formula which takes account of these bifurcations, and we examined its semiclassical 
accuracy. 

In the second part of this paper we deformed the ellipse into an oval which is non-integrable. It 
was shown that in the quasi-integrable regime it is not always sufficient to treat the break-up of tori 



by the general formulas of [18, 2C, If a torus in the integrable system is close to a bifurcation, 
then uniform formulas for general bifurcations are needed for describing its break-up. In the numerical 
section it was shown that with the uniform formulas for the break-up of tori and for the bifurcations 
the same semiclassical accuracy can be achieved as in the unperturbed integrable system. 

If the billiard system is deformed more, then the stable periodic orbits of the billiard undergo 
further bifurcations. A semiclassical description of a system in the truly mixed regime, i. e. not in the 
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Figure 9: The semiclassical errors D QM {x) — D QTI (x) (dotted) and D QM (x) — D UA (x) (full) for the 
quasi-torus approximation with interference term and for the uniform approximation, respectively. 
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Figure 10: The semiclassical errors D QM (x) — D UA (x) (dotted) and D QM (x) — D CUA {x) (full) for the 
uniform approximation with and without the contribution of a pair of complex orbits, respectively. 
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near-integrable regime, requires the use of uniform approximations for all generic bifurcations ^9|. In 
case that the system has discrete symmetries there are further bifurcations which are specific for the 
considered symmetry. The application of uniform approximations for these bifurcations requires not 
only the determination of all real periodic orbits up to some length, but also of all complex periodic 
orbits which are close to becoming real. Additional complications can arise when a periodic orbit 
undergoes several subsequent bifurcations which cannot be considered isolated. This shows, as is 
well-known, that semiclassical trace formulas for systems with mixed phase space are definitely more 
complicated as for integrable or chaotic systems. 
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